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Abstract 

A monochromatic point source, embedded in a three-dimensional disordered 

medium, is considered. The resuhing intensity pattern exhibits a new type of 

long-range correlations. The range of these correlations is infinite and their 

magnitude, normalized to the average intensity, is of order l/koi, where /cq 

and i are the wave number and the mean free path respectively. 
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A wave propagating in a random medium undergoes multiple scattering and produces 
a complicated, irregular intensity pattern. Such a pattern is described in statistical terms 
and one of its important characteristics is the correlation function (/(ri)/(r2)), where /(r) 
is the wave intensity at point r and the angular brackets designate the ensemble average. 
This correlation function has been studied in recent years, both theoretically and experi- 
mentally, and various types of correlations has been identified^. There are large short-range 
correlations, due to interference among the waves arriving, via all possible scattering se- 
quences, to a neighborhood of a given point; this interference is responsible for the large 
spatial fluctuations in intensity (speckles). There exist also weak long-range correlations, 
due to diffusion which propagates the locally large fluctuations to distant regions in space. 
The purpose of this note is to identify and study a new type of long-range correlation which, 
under appropriate circumstances, can dominate the "standard" long-range correlations^. 

I consider a monochromatic point source embedded in an infinite three-dimensional ran- 
dom medium. The position of the source is tq and its frequency is u. Assuming a scalar 
wave, one has the following equation for the field (the Green's function) at point r: 

{V^ + kl[l + fxir)]+iri}G^if,fo)=5if-fo) . (1) 

Here fi{r) is the fractional fluctuation of the dielectric constant, r^ is a positive infinitesimal 
and ko = w/c, where c is the speed of propagation in the average medium. I assume that 
/x(r) obeys white-noise Gaussian statistics, i.e. 

(/x(f)/x(f')) = u5{f- r*) , (2) 

where the constant u describes the strength of the disorder. 

For weak disorder, the average field, the average intensity and various correlation func- 
tions can be computed by the diagram technique^. The average field, (G^(r*, tq)), decays 
exponentially away from the source: 



{G^{f, fo)) = -^ — ^ exp 
47r I r — ro I 



(^^0-^) I r-fo 



(3) 



where d, = Air/ukQ is the mean free path and koi ^ 1. 

The intensity at point r is defined as /a;(^) =| G^{f,fQ) p and its average value is given 
by the diagram in Fig. 1. Intensity propagates from the source to a distant observation 
point, such that | r — ro |^ £, by a diffusion process. This is represented by a diffusion 
ladder (the shaded box), T(ri,r2) = 3/i^ | ^i — "^2 I- The vertices, connecting the external 
points to the ladder, are short-range objects which can be replaced by the number 

rfVi(G.(fo,n))(G:(fo,fi)) = ^ , (4) 

so that the average intensity at point r is 

(Url) = if] nf,ro) = ., ,,|l -. I . (5) 
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Let us now consider the correlation function 

{AUr-)AUr + Af)) ^ 

where A/(^(r) = /(^(r) — (/(^(r)) is the deviation from the average value. The points r and 
r + Ar are assumed to be far from the source, i.e., many mean free paths away. Since only 
long-range correlations are considered here, Ar ^ £ is assumed. 

The leading contribution to the long-range correlations is represented by the diagram in 
Fig. 2. The diagram contains two pairs of Green's functions. One pair propagates, via a 
diffusion ladder, to point f; the other - to point f+Af. The pairs are connected by a dashed 
line, which makes this diagram belong to the set of diagrams contributing to C(Ar). The 
physical interpretation of this diagram is straightforward: the wave, emanating from the 
source, gets scattered at some point fi, close to the source. The secondary wave, emerging 
from fi, propagates by diffusion to two distant points, r and f+ Af. In this way intensity 
correlation at the two points is established. 

In order to evaluate the diagram in Fig. 2 one has to compute the vertex, depicted 
separately in Fig. 3. The vertex is a short-range object and, thus, can be replaced by a 
point, located at rg, from which the two diffusion ladders emerge. The number V, assigned 
to the point is: 






x(G^(r-'3-fi))(G^(fi-fo))(G:(f3-fo)) . (7) 

Integration over r2 gives 

d'MG^if, - fo))(G:(f2 - n)) = ^U\ ro -n\). (8) 



47r 
where 



sin /coX f x\ /n^ 

Integration over f^ gives an identical contribution, so that 

V = y{'L) Jd'n{GAn-ro)){G:{n-fo))a\ fo-n I) = 

The integral is dominated by the region of small p, i.e. p ~ k^^ <^ £, so that the exponential 
can be replaced by 1, which leads to V = i/32iTkQ. The diagram in Fig. 2 is, thus, equal to 

T(fo,f)T(fo,f+Af) — (11) 



where the factor f ^ j accounts for the two vertices connecting the diffusion ladders to the 
observation points r and f+ Ar. Recalling that (/^^(r)) = (£/47r)^T(r, ro), and similarly 
for {Iuj{r + Ar)), one can rewrite Eq. (11) as (7r/2A;o^)(/aj(r))(/tj(r + Ar)). Finally, the 
diagram in Fig. 2 should be assigned a combinatorial factor 2, since one could connect the 
two external Green's function by a dashed line, instead of connecting the two internal ones 
as in Fig. 2. (Only lines with arrows in opposite directions, corresponding to a pair of G, G* 
need to be connected.) Thus, the contribution of the diagram in Fig. 2 to the normalized 
correlation function [Eq. (6)] is: 

C„(A.) ^ ^ (12) 

where the notation Co is used to distinguish this term from the other types of correlations, 
usually designated^ as Ci, C2 and C3. For the geometry considered, i.e. point source in 



an infinite three-dimensional medium, the Co-term is much larger than the other types of 
long-range correlations. Moreover, it does not decay in space (infinite-range correlation). 

It is worth to note that the Co-term resembles the C2-term^'^, in the sense that both 
terms descibe correlation between two distant points, for intensity propagating from a single 
source. However, Co is not contained in C2. Indeed, C2 is desribed by a 4-ladder diagram 
and, thus, cannot contain the diagram in Fig. 2. Moreover, in computing and interpreting 
the C2-correlation it is essential that all 4 ladders are "at work" and can be treated in the 
diffusion approximation^'^. For a point source in three dimensions the C2-term is of order 
(/co^)"^, which is much smaller than the Co-term in Eq. (12). 

So far Ar ^ d. was assumed. The case Ar = 0, however, is also of interest. For this 
case the diagram in Fig. 2 gives a contribution to the second moment of intensity. For 
the normalized intensity, /(r) = I{r)/{I{r)), this contribution is 2Co(0) = 271 /k^d, (the 
extra factor 2, as compared to Eq. (12), appears because for Ar = one can pair also 
Green's functions emerging from point r). This contribution represents a small correction 
to the Rayleigh value {P) = 2. The knowledge of the correction enables one to compute 
deviations of the intensity distribution P{I) from the Rayleigh function Po(-^) = exp(— /). 
The calculation is practically identical to the one performed by Mirlin et al.^, with one 
important difference: Ref. 3 considered quasi-one dimensional geometry, where the correction 
was quite different from the present value 27r/koi. For the present, three-dimensional case 
the result is: 

P(/) ^ exp (-/ + ^P) , / « koi . (13) 

The asymptotic tail of the distribution, for / ^ k^i, can be obtained by the method of 
optimal fluctuation*^. For / ~ k^i the tail should smoothly match expression (13). A 
detailed calculation of P{I) will be presented in a separate publication. 

One can easily extend the calculation to include correlations at different frequencies. The 
object to consider is the same as in Eq. (6), but with Iuj+Auj{r + Ar) instead of /(^(r + Ar). 
The only difference in the calculation is that, in Eq. (7) for V, one has to substitute uj + Au 



for oo in the last three Green's functions. This leads to a replacement 00^00 + Auj in one of 
the Green's function and in one of the f-factors in Eq. (10). For Auj ^ uj this replacement 
leads only to inessential correction to the constant value in Eq. (12). 

In conclusion: 
(i) It has been shown that the intensity pattern, produced by a wave propagating in a 
random medium, can exhibit a new type of correlation, designated as Cq. These correlation 
is of order 1/koi and has infinite range in space and frequency. Unlike the other types of 
long-range correlations, C2 and C3, the Co-term is dominated by the configuration of the 
disorder near the source. 

(ii) Being sensitive to the short- distance properties of the disorder, the Co-term cannot be 
universal, i.e. its magnitude should depend on the specific type of the disorder. For the 
white-noise Gaussian disorder, considered in this note, Cq is determined by the parame- 
ter koi. However, for Gaussian disorder with a finite correlation length d, or for discrete 
impurities of size d, the Co-term should depend on d. 

(iii) Let us emphasize that this paper considers a propagating wave, in an open system. A 
somewhat different problem pertains to the statistics of eigenmodes, and eigenfrequencies, 
in a closed disordered system. Non-universal terms, of the same origin as these considered 
here, exist also in that case, as has been recently discussed by Mirlin^. 
(iv) The Co term is also present in lower dimensionalities. The quasi-one dimensional ge- 
ometry (a tube of length L and cross section A) is of a particular interest, since it is often 
used in experiments^. The C2 and Cs-terms in this geometry are of order g^^ and g^^, 
respectively, where g ~ kQAi/L. Thus, for a point source inside the tube, the Co term will 
dominate as long as koi < g, i.e. L < Ako. This implies that our earlier calculation^, which 
left out the Co-term, is valid only for L > Ak^. 
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the present work. Useful conversations with A. Genack, D. Khmelnitskii and A. Lagendijk 
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